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ABSTRACT 

Using the SARA 0.9 meter telescope at Kitt Peak National Observatory (KPNO) we obtained R 
images during the summer of 2010 and V images during the summer of 2012 of the globular c luster 
M14. These images were analyzed using the DAOPHOT/ALLFRAME package of |Stetson| ([19871 



1994 1 to create a preliminary color-magnitude diagram (CMD) of M14. We have positively identified 
the positions of 64 of the 112 known RR Lyrae stars on the CMD. We find considerable spread in 
both the magnitudes and color of the RR Lyrae stars indicating a significant amount of differential 
reddening on relatively small scales. This is also seen in the giant branch which shows up to a ^0.15 
magnitude spread in V — R color. We also found the median magnitude of the RR Lyrae stars and 
hence the horizontal branch to be V—17.2. 

Subject headings: clusters: globular, stars: variables: general-Galaxy: individual: M14 



1. INTRODUCTION 

Globular clusters are relatively old and concentrated 
star clusters found in the galactic halo. Because the stars 
in a given cluster are presumed to have formed from the 
same cloud of gas and dust they most likely have simi- 
lar composition, age, and distance. This makes globular 
clusters ideal laboratories for studying the evolution of a 
single population of stars. Depending on their metallic- 
ity, age, and mass globular clusters can have large num- 
bers of RR Lyrae stars. RR Lyrae stars are pulsating 
variable stars located on the horizontal branch with ap- 
proximate absolute magnitudes of U=0.5, making them 
ideal standard candles. They come in two general types, 
RRO variables which have longer periods and larger am- 
plitudes than RR1 variables which have a more sinusoidal 
light curve. Though they have very similar V absolute 
magnitudes RR1 variables are slightly bluer than their 
RRO counterparts, differing in color by roughly 0.1 in 
(V — R). By using Fourier decomposition of the variable 
light curves it is possible to compare the amplitudes and 
phases of the light curves to hydrodynamic stellar evo- 
lution models. From these comparisons estimates of the 
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mctallicity, h elium abundances, and ages of these sta rs 
can be made ( Simon k, Teays|[l982 Simon "et~al]p93| . 

Early studies of variable stars m globular clusters were 
made using photographic plates, and despite the errors 
involved in their use many variables were found, the most 
comm on being of RR Lyrae type (Wehlau & Froelich 
1994 1. However, even to this day the cores of many glob- 



ular clusters are yet to be examined using modern tech- 
niques to search for variables. With the use of CCDs 
and image subtraction software it is possible to discover 



heretofore unknown variable stars ( Alard & Lupton 1998 



Alard 20001). In many cases the number of known vari 



abies in a cluster may double using these new techniques. 
Having a larger sample of variables can allow for a more 
detailed analysis of a cluster and a better determination 
of its physical parameters. 

Previous work by our research group has focused on 
several globular clusters, M14 (NGC 6402), NGC 4833 
NGC 6584, and more recently M107 (NGC 6171) flCon- 
roy et al.||20l2l |Darragh fc Murphy1[20T2| |Toddy et al 



2012} |McCombs et al.||2012|>. IrTtTii s study we revisit 
M14. Recently |Conroy et al. (2012), using the image 
subtraction method, found 59 previously unknown RR 
Lyrae variables in M14, more than doubling the total 
and raising the number to 112. Because of this relatively 
large population of RR Lyrae variables M14 became an 
ideal case for further study. 

M14 is relatively close to the galactic center (/ = 21.3°, 
b = 14.8°). Because of this proximity M14 is signifi- 
cantly reddened due to interstellar dust and has a listed 
E(B — V) of 0.6 implying nearly 2 magnitudes of ex- 
tinction in V (Harris 1996). It is also likely that this 
extinction is not uniform but instead patchy across the 
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Fig. 1. — 100 micron emission in the vicinity of the globular 
cluster M14 (Schlegel, Finkbeiner, & Davis 1998). This is a repro- 
cessed composite of the COBE/DIRBE and IRAS/ ISSA maps, 
with the point sources removed. The contours cover the range 
from 9.01 (blue) to 22.31 MJy/sr (red) in intervals of 2.22 MJy/sr. 
The square near the center of this image represents our field of 
view with SARA KPNO. 

cluster. Figure 1 shows a 5 x 5° region at 100 microns 
centered on M14. It is clear that there is significant dif- 
ferential reddening across the region. This differential 
reddening can prevent an accurate determination of the 
distance, age, and metallicity of the cluster. It is also a 
concern when constructing a CMD of the cluster. In spite 
of this difficulty we have used data taken with the South- 
eastern Association for Research in Astronomy (SARA) 
telescopes to create a CMD and identify the location of 
RR Lyrae stars on it. From these CMD positions we 
are then able to verify if they are consistent with cluster 
membership. The CMD also allows us to estimate the 
amount of differential reddening in the vicinity of M14. 

2. OBSERVATIONS & REDUCTION 

i ?-band image frames u sed in this study were published 
by Conroy et al. (2012). Those observations were ob- 
tained in May and June of 2010 using the KPNO SARA 
0.9 meter telescope with an Apogee Alta U42 CCD with 
a 2048x2048 Kodak e2V CC42-40 chip, a gain of 1.2 elec- 
trons per count, RMS noise of 6.3 electrons, and cooled 
to a temperature of approximately — 30°C. lxl bin- 
ning was used, resulting in a scale of 0.42"/px and a 
13.6x13.6' field of view. Typical seeing was 2.2" and 
ranged from 1.5-3.0". Exposure times were set at a con- 
stant 60 seconds to avoid overexposure of the core or 
bright giants during best possible seeing conditions. V- 
band images were obtained during the summer of 2012 
also using the SARA 0.9 meter telescope on the nights 
of 14 May and 7, 8, 9 June 2012. The same Apogee 
Alta U42 CCD was used with an exposure time of 100 
seconds, however we used 2x2 binning. This resulted 
in a scale of 0.825"/px and a 13.6x13.6' field of view. 
We used standard data reduction procedures using the 
Maxim DL image processing software. Because hot or 



dead pixels as well as cosmic rays can produce a large 
number of false positives when running the image sub- 
traction software ISIS, we also applied a bad pixel map 
and ran Maxim DL's hot and dead pixel kernel filters. 
We then applied the CRAVERAGE cosmic ray removal util- 
ity in IRAF to each frame. 

3. ANALYSIS 

3.1. Image Subtraction 

ISIS v2.2 was used to search for, identify, and deter- 
mine the positions of the variables (Allard & Lupton 
1998; Allard 2000). The package contains seven c-shell 
scripts and three parameter files that were used in our 
analysis. ISIS uses pattern recognition to align stars from 
image to image and search for changes in brightness of 
every star it finds. ISIS creates a reference image from a 
set of the highest quality images so that all other images 
are compared/subtracted from this. ISIS also takes into 
account differences in seeing. To do this ISIS convolves 
the reference image so that the point-spread function of 
the stars is nearly identical to that of the image it is be- 
ing subtracted from. Note that ISIS uses relative flux 
not magnitudes. The ISIS parameters and methods that 
we used are described in m ore detail in the analysis of 
M14 bylConroy et all (|2012h. 



3.2. Photometry 

While ISIS very useful for searching for variables it 
does not allow for absolute photometry of the vari- 
ables. The relative flux output from ISIS is not ade- 
quate since the magnitude scale is not linear. The rel- 
ative fluxes are useful for determining period and vari- 
able type but of little use in determining most photo- 
metric properties of variables. Therefore it is necessary 
to obtain instrumental magnitudes of each of the vari- 
ables found in addition to other stars in the cluster so 
that a CMD can be created and cluster membership con- 
firmed. Another added difficulty are the crowded fields 
typically found in globular clusters, particularly their 
cores. To deal with the crowded nature of globular clus- 
ters it has become standar d procedure to use S tetson's 



DAOPHOT / ALLFRAME ( |Stetson||i987 |1994[ ) to cre- 
ate CMDs. The combination of these packages results in 
magnitudes for both V a nd R that we then could match 
using a script written by Tabur (2007). 



3.3. Magnitude Transformation 

ALLFRAME gives instrumental magnitudes of each 
star in the images with an arbitrary zero point of 25 and 
does not take into account exposure time, extinction, or 
color. We transformed the instrumental magnitude Vi nst 
to true apparent magnitude V app by using a set of stan- 
dard stars within the cluster as identified and calibrated 



by Stetson ( 2000 1 . These standard stars covered a wide 
range of magnitudes that resulted in a linear relationship 
between our instrumental magnitude and Stetson's stan- 
dards. The next step was to transform the instrumental 
(V — R) to true values. This was done using transfor- 
mation coefficients determined earlier in 2012 by one of 
the authors. It is also worth mentioning that because of 
the differing binning between the V and R images the R 
images were converted from lxl to 2x2 so that proper 
alignment and colors could be determined. This resulted 
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Fig. 2. — Color Magnitude Diagram of M14 from two V and R. 464 images were stacked together to produce the V image and 698 
images were stacked to produce the R image. The blue points are the RR1 stars and the red are the RRO. Much of the scatter in the RR 
lyre stars is due the result of stacking images to improve signal to noise and due to differential reddening. The diagonal lines indicate the 
reddening tracks. 



in counts that were 4 times less than would be expected. 
The factor of 4 difference occurs when Maxim DL does 
pixel operations of this nature. Maxim DL interpolates 
over the pixels and point spread functions then assigns a 
newly binned pixel value that is equal to the median of 
the original pixels. In reality the new pixel value should 
be proportional to the increase in area of the new 2x2 
rebinned pixel, i.e., a factor of 4. To correct for this we 
multiplied each pixel in the R frame by 4. 

4. RESULTS 

The color magnitude diagram of M14 was created using 
V and R band KPNO data. The V and R average images 
were created by stacking images. The V stacked image 
was an average of 464 images and the R stacked image 
was created using an average of 698 images. Because of 
this stacking the magnitudes of our RR Lyrae variables 
are not average values but dependent upon when they 
were observed. This results in some scatter of these vari- 



able stars in the CMD. Stars in the V and R s tacked 
images were matched using a script (Tabur 2007). The 



script creates a transform between the coordinates of the 
stars detected in the V and R stacked frames, and only 
retains those stars that are detected in both colors. Thus 
any stars that did not have matches in each color were 
eliminated by the script. This process resulted in 4721 
star matches in both V and R that were then be plotted 
in the CMD. 

ISIS and ALLFRAME both give positions of variables 
and all stars found, respectively. But because ISIS used 
only a few select good seeing frames that were combined 
into a reference image for subtraction and ALLFRAME 
used a much larger number of stacked frames, the re- 
sulting positions of the variables found by ISIS can vary 
somewhat from ALLFRAME. To ensure that we were 
properly pulling our variables out of ALLFRAME we 
only plotted variables found by ISIS that were within 0.5" 
of the ALLFRAME position. That lowered the scatter of 
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the variables on the CMD and eliminated any misiden- 
tifications. In total we were able to identify and plot 33 
RRO and 31 RR1 variables onto the CMD. 

The CMD for M14 is shown in Figure 2. The RR Lyrae 
stars are shown by the colored dots, with the RRO's being 
red and the RRl's shown in blue. The diagonal lines indi- 
cate the direction of reddening. It is apparent that there 
is significant scatter in the diagram. Most of scatter is 
due to differential reddening. Compared to other clusters 
observed by this group the scattered found here is much 
greater than say for NGC 6584 whic h is much less red- 
dened than M14 (|Toddy et aT]|2012[). As expected the 



largest scatter is found among the RRO variables since 
they have the largest amplitudes and greatest change in 
color from minimum to maximum phase. The scatter is 
somewhat less for the RR1 variables due to their smaller 
amplitudes. Given their shorter periods of roughly 8 
hours we were more likely to be close to their time- 
averaged magnitudes when stacking images to produce 
the CMD. The two RR1 variables found well above the 
horizontal branch are likely blended with redder stars. 
Several of the RRO variables to the upper left of the in- 
stability strip were probably those that happened to be 
imaged each night near maximum. At peak brightness 
RRO variables are much bluer than average, thus this 
scatter can mimic reddening. We found the median V 
magnitudes and (V — R) color of the RRO variables to be 
17.26 and 0.91, respectively. For the RR1 variables we 
found these quantities to be 17.10 and 0.82, respectively. 
The color difference is typical given that RR1 variables 
are warmer than their RRO counterparts. Using all 64 of 
these RR Lyrae variables we find a horizontal branch V 
magnitude of approximately 17.17, similar to but slight ly 
brighter than the cataloged value of 17.3 ( |Harris|[l996[ ). 

By examining the giant branch of the cluster it is possi- 
ble to estimate the degree of differential reddening. The 



A(V — R) is roughly 0.1 to 0.15 from one side of the giant 
branch to the other. This is greater than the expected 
photometric errors. Given this spread of 0.15 magnitudes 
in color the expected AV could be 0.5 to 0.75 magni- 
tudes. This is close to what is seen in the V variation 
among the RR Lyrae variables, indicating significant dif- 
ferential reddening as would have been expected from the 
100 micron image shown Figure 1. 

5. CONCLUSIONS 

Due to time constraints, we have only presented a pre- 
liminary analysis of the the CMD and location of the RR 
Lyrae variables in M14. From the 64 positively identified 
variables on the CMD we estimate a horizontal branch V 
magnitude of 17.2. Our analysis has shown there is sig- 
nificant scatter in the RR Lyrae variables, much of it due 
to differential reddening that must be addressed in our 
future investigations in order to determine a more accu- 
rate distance, age and metallicity of the cluster M14. To 
accomplish this our next task will be to determine the 
time-weighted V and R magnitudes for the RR Lyrae 
variables via full Fourier decomposition. Once this is 
completed we will be better able to determine physical 
parameters of the RR Lyrae variables and those of M14. 
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